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Abstract 

We present quantum mechanical close-coupling calculations of collisions between two hydrogen 
molecules over a wide range of energies, extending from the ultracold limit to the super-thermal 
region. The two most recently published potential energy surfaces for the H2-H2 complex, the 
so-called DJ (Diep and Johnson, 2000) and BMKP (Boothroyd et al., 2002) surfaces, are quan- 
titatively evaluated and compared through the investigation of rotational transitions in H2+H2 
collisions within rigid rotor approximation. The BMKP surface is expected to be an improve- 
ment, approaching chemical accuracy, over all conformations of the potential energy surface com- 
pared to previous calculations of H2-H2 interaction. We found significant differences in rotational 
excitation/de-excitation cross sections computed on the two surfaces in collisions between two para- 
H2 molecules. The discrepancy persists over a large range of energies from the ultracold regime 
to thermal energies and occurs for several low- lying initial rotational levels. Good agreement is 
found with experiment (Mate et al., 2005) for the lowest rotational excitation process, but only 
with the use of the DJ potential. Rate coefficients computed with the BMKP potential are an 
order of magnitude smaller. 



I. INTRODUCTION 



Collisions involving two hydrogen molecules are of great interest for three main reasons. 
First, the H2+H2 collision system is a prototype for chemical dynamics studies and can 
be used as a testing ground for scattering theory of inelastic (non-reactive) diatom-diatom 
collisions involving a weak interaction potential. Second, H 2 is the most abundant molec- 
ular species in the universe. The rotational and vibrational transitions in H2 induced by 
collisions with its twin are of practical importance in models of astrophysical environments 
where the physical conditions may not be accessible to terrestrial experiments. Examples in- 
clude low densities characteristic of giant molecular clouds in the interstellar medium where 
star formation occurs and H2 may act as a coolant 1 . Heating of the interstellar cloud by 
strong shock waves induces rotational and vibrational excitation of the H 2 molecules leading 
to collision-induced dissociation to two free H atoma 2 - and photodissociation regions where 
H 2 is exposed to strong UV stellar radiation 3 . Third, with recent experimental advances 
in the cooling and trapping of molecules^^ to form molecular Bose-Einstein condensates, 
collisional studies of the H 2 +H 2 system can serve as a model to provide new insights into 
the behavior of diatom-diatom collisions at ultracold temperatures including investigations 
of Feshbach resonances, predissociation in van der Waals complexes, determination of com- 
plex scattering lengths, testing of effective range theory and Wigner threshold laws, and 
quasiresonant vibration-rotation energy transfe r^^ 1 ^ 11 ' 12 in molecular collisions. 

Since (H 2 )2 is the simplest 4-identical-particle closed shell system and there is a con- 
tinuing demand for accurate collisional data for modeling astrophysical, atmospheric and 
combustion processes there have been considerable theoretica l 13ll4ll5ll6ll7ll8ll9l2 ° l21122123124125 
and experimenta l 2 ^ 2 ? 12 ^ 2 ? 13 ' 1 ' 131 studies performed on the H2+H2 system. The most recent 
experiment was performed by Mate and co-workers 31 . Using the technique of Raman spec- 
troscopy with supersonic expansions of para-H 2 , they measured the rate coefficient, &oo->2o> 
for H 2 (ji = 0) + H 2 (j2 = 0) — > H 2 (ji = 2) + H 2 (j 2 = 0) collisions in the temperature range 
of 2 to 110 K. With their experimental methodology, without the loss of generality, the 
reduction of the measured /coo^2o rate coefficient to the corresponding cross section, 000^20 
in the collision energy range of 360 to 600 cm -1 was made possible. To complement the ex- 
perimental measurement, Mate et alQ also performed coupled channel quantum scattering 
calculations to determine the cxoo^2o excitation cross section as a function of the collision 



energy. Good agreement was found between the experimentally-derived cross sections and 
theoretical results obtained using the rigid rotor potential energy surface (PES) developed 
by Diep and Johnson (DJ)^2 

The most recent theoretical study of rotational excitation in H 2 +H 2 collisions was car- 
ried out by Gatti and co-workers^. They employed the wave-packet propagation method 
in conjunction with the multiconfiguration time-dependent Hartree algorithm to compute 
rotational excitation cross sections for collision energies up to 1.2 eV by a flux analysis of 
the interaction of the wave-packet with a complex absorbing potential. Gatti et al. com- 
pared their results with the wave packet calculations of Lin and Guo^i who employed the 
coupled-states (CS) approximation which neglects Coriolis coupling. They found that the 
CS approximation can lead to rather reliable predictions provided the calculations are per- 
formed at low collision energies and low rotational excitations. However, without the aid 
of accurate molecular potential curves, it is difficult to determine in which energy regime 
the Coriolis coupling takes effect for a specific collision system. Both wave-packet calcula- 
tions employed the PES constructed in 2002 by Boothroyd, Martin, Keogh, and Peterson^ 
(BMKP). While previous calculations^ have indicated that the CS approximation can give 
accurate results for rotational excitation in H2 at low energies, these calculations have used 
rather simple potential functions for the H4 system. 

The purpose of this study is to perform accurate close-coupling calculations of rotational 
transitions in H2+H2 collisions within the ground vibrational state. Due to the relatively 
large vibrational level spacings of the H 2 molecule and the weak H 2 -H 2 interaction poten- 
tial, the rigid rotor approximation is expected to hold well for the present system for pure 
rotational energy transfer. This is also confirmed by the wave packet calculations of Lin 
and Guo^. Here, we perform close-coupling calculations of rotational excitation on both 
the DJ and the BMKP PESs and compare our results with the wave-packet results of Lin 
and Guo and Gatti et al. So far, the accuracy of the BMKP PES has not been established 
by comparing results from accurate quantum calculations with experimental results. Such 
a comparison is provided here for both the DJ and BMKP PESs and we show that the rigid 
rotor DJ potential yield results that are generally in better agreement with experiment. 

We also investigate the behavior of elastic and rotationally inelastic collisions in diatom- 
diatom collisions at ultracold energies by taking the H 2 +H 2 system as an illustrative example. 
While ultracold rotational and vibrational energy transfer in atom-diatom collisions have 



extensively been reporte d^^ 1( - ) i 11 i 12 i^ in the last several years, such calculations are yet to 
be performed on molecule-molecule collisions. ForreyJ^ has performed limited calculations 
on H 2 +H 2 collisions in the ultracold regime using the semi-empirical potential of Zarur and 
Rabita^ while Avdeenkov and Boh n 37 i 38 i 39 reported spin-exchange collisions in O2+O2 and 
OH+OH/OD+OD systems. Here, we provide a detailed investigation of rotational energy 
transfer in H 2 +H 2 collisions on the DJ and BMKP PESs at ultracold energies and report 
complex scattering lengths for collisions involving different initial rotational levels of the two 
molecules. 

The paper is organized as follows: A brief description of the theoretical methodology is 
given in section |H] and results are presented in section IIHI Section IIVI provides summary 
and conclusions. Atomic units are used throughout, unless otherwise noted: i.e., h — e — 
m e = a Q = 1 a.u., while 1 hartree = 27.2116 eV = 627.51 kcal/mol. 



II. THEORY 



Calculations of state-to-state rotational transition cross sections and rate coefficients can 
provide an important test of the reliability of the potential energy surfaces describing the 
interaction of two H 2 molecules when compared to available experiments. To compute the 
scattering amplitudes and hence the cross sections, we use well established quantum mechan- 
ical close-coupling method o 14 i 17 i 40 i 41 . Here we only summarize the essence of the theory. To 
describe the scattering process, we solve the time-independent Schrodinger equation within 
rigid rotor approximation for the collision of two H 2 molecules in the center of mass frame 
given by 

(f(R) + J2 hrotih) + V{R, r u f 2 ) - E^j V{R, h, h) = 0, (1) 

where T(R)=—-^S7 2 R is the kinetic energy operator and h rot {ri) = -^-pi is the diatom rota- 
tional kinetic energy operator; fi and fii are the reduced masses of the H 2 collision pair and 
a isolated H 2 molecule, respectively. The internuclear distance between the two H atoms is 
denoted by r^, and R is the distance between the center of mass of the diatoms; f 1 and f 2 are 
the orientation angles of the rotors 1 and 2, respectively. The term V(R, r\,r2) represents 
the H 2 -H 2 interaction potential. 



The rotational part of the Hamiltonian operator satisfies the eigenvalue equation 

(hrotin) - Bijiiji + 1)) Y jim . (ti) = (2) 

where Bi are the rotational constants of the rigid rotors. 

The rotational angular momenta of the two molecules j\ and j 2 are coupled to form 
ji2, which is subsequently coupled to the orbital momentum I to yield the total angular 
momentum J. The angular wave function in the total angular momentum representation is 
given by 



Yj 2 m2 {h)Yi mi {R) (3) 

and under spatial inversion 

P#2"(M,*a) = {-l) h+h+h ^ J f {R,h,r 2 ) (4) 

where a = jij 2 ju, and mi, m 2 , m\ 2 and m\ are the projections of j\, j 2 , j\ 2 and I, respec- 
tively, onto the space-fixed z-axis. The symbol (jij 2 mim 2 \JM) denotes a Clebsch-Gordon 
coefficient. 

Using the basis set expansion method, we expand our total wave function ansatz as 

m,h,r 2 ) = F™W"lfrh>h) (5) 

JMal 

Substituting eqn.© into eqn.((TJ), we arrive at a system of close-coupling equations 

" + 2fiEk ) Fil{R) = 2 ^ F ^(R)(€f\v\4> J Jf'), (6) 

^ ' a'l' 

where the quantities in (• • • ) are the matrix elements of the interaction potential and E^ = 
E — Biji(ji + 1) — B 2 j 2 (j 2 + 1) is the kinetic energy of the relative motion for a given value 
of the total energy E. The solution of the coupled equations and asymptotic analysis of the 
radial wave functions yield the scattering S-matrix from which cross sections for state-to- 
state rotational transitions from an initial level specified by quantum numbers jij 2 to final 
levels j[j' 2 are given by 



where v = jij 2 ji 2 l- 

The above expression for cross section assumes that the two diatomic molecules are 
distinguishable. However, for H 2 -H 2 collisions the target and projectile molecules are in- 
distinguishable and one must take into account the symmetry of the wave function under 
exchange. Thus, properly symmetrized total angular momentum wave functions 17 

±(-l)^+^+V/ 1 ? 2 ^(4ri,r 2 )] (8) 

need to be employed in which the index pair jij 2 is restricted to ji > j 2 to obtain a 
linearly independent set. Using the symmetrized angular wave functions, one obtains coupled 
equations similar to eqn. (jHJ) which yield scattering cross sections^ 7 - 

7r(l + )(T + 8 Alt' ) T n 

Rate coefficients for state-to-state rotational transitions are obtained by averaging the 
appropriate cross sections over a Boltzmann distribution of relative speeds of the projectile 
molecule at a given temperature T: 

poo 

k jlh ^ m (T) = G / dE k a hj ^{E k )E h e^\ (10) 
Jo 



where the constant G = yj^f3 2 and (3 = (ksT) 1 with ks being the Boltzmann constant. 
The total quenching rate coefficient can be calculated from 

k hh ( T ) = Yl ( T ) • ( 1 1 ) 

III. RESULTS 

We have carried out close-coupling calculations for collisions of H 2 with H 2 using the 
BMKP and DJ PESs. The rigid rotor target and projectile energy levels were calculated 
using a rotational constant of B = 60.853 cm -1 for the H 2 molecule. To solve the coupled 
radial equations (JBJ), we used the hybrid modified log-derivative-Airy propagatoi^ in the 
general purpose non- reactive scattering code MOLSCAT 44 . The log-derivative matrhj 4 ^ is 
propagated to large intermolecular separations where the numerical results are matched to 



the known asymptotic solutions to extract the physical scattering matrix. This procedure is 
carried out for each partial wave until a converged cross section is reached. We have checked 
that the results are converged with respect to the number of partial waves as well as the 
matching radius for all channels included in the calculations. 

In addition to the partial wave convergence, based on the DJ PES, we have checked that 
the results are converged with respect to various parameters that enter into the close-coupling 
calculations. These include the number of quadrature points used for angular integration, 
the number of terms in the angular expansion of interaction potential, and the asymptotic 
matching radius for radial integration. In the ultracold regime, we used a matching radius 
of R Q = 200-300 a.u. to obtain converged values of elastic and inelastic cross sections while 
for the subthermal energy region, a matching radius of 50 a.u. was sufficient to yield results 
of comparable accuracy. Similarly, we used 10 quadrature points each for integration along 
angular coordinates 9±, 9 2 , and 0i 2 . 

Finally, two different basis sets (22-state: = 00, 20, 22, • • •, 44 and 50-state: jij 2 
= 00, 20, 22,- • •, 66) were also employed to further test the convergence of our results. For 
E < 1.0 eV, the two basis sets yield results within 1%, and at E — 1.0 eV, a similar degree 
of accuracy was obtained for the dominant transitions. For weaker transitions such as 00 
— > 44 rotational excitation, the small and large basis sets gave a o^44 = 3.01 xl0~ 18 and 
3.20 xl0~ 18 cm 2 , respectively, at 1.0 eV. Note that the cross sections for dominant transitions 
are two orders of magnitude larger than the weaker ones. Even at E — 2.6 eV, the cross 
sections obtained from both basis sets for the dominant transitions have similar convergence 
properties as for E — 1.0 eV. However, for 00— >44, we found that the larger basis set is 
preferred since there is a significant difference between the cross sections (i.e., 000^44 = 
4.31 xlO" 17 cm 2 and 6.30 xl0~ 17 cm 2 , respectively, from the small and large basis sets, at 
2.6 eV). Since our focus is on the low-energy region where the two basis sets yield similar 
results, the smaller basis set is adopted throughout the calculations. 

Fig. 1 shows the comparison between theoretical and experimental rate coefficients for the 
00—^20 transition in the temperature range between 50 K and 300 K. Both experimental and 
theoretical rate coefficients indicate a precipitous drop for temperatures lower than 100 K. 
Unexpectedly, we find that the 00— >20 excitation rate coefficient computed with the BMKP 
PES is about an order of magnitude smaller than that calculated with the PES of D J and the 
experiment, though both BMKP and DJ display the same trend. Only results from the DJ 



PES agree with the experimental data of Mate et alQ Good agreement between experiment 
and theory based on the DJ PES was also shown by Mate et alQ. Further, the theoretical 
results obtained by Flower— using the older PES of Schwenke^ are also seen to be in good 
agreement with experiment. The discrepancy with the BMKP results may directly be traced 
to the weaker anisotropy of the BMKP PES responsible for the 00— >20 transition. In the 
calculations the angular dependence of the interaction potential is represented as^ 

V(R, 9 U 2 , 12 ) = V hM {R)G hM {9 x , 9 2 , 12 ), (12) 

where Vi u i 2t i(R) are radial expansion coefficients and Gi lt i 2 j(9i, 9 2 , (fin) are bispherical har- 
monics. In Fig. 2 we compare the spherically symmetric (Vooo) and the leading anisotropic 
terms, V022 = V202 and V224, in the angular expansion of the BMKP and DJ interaction 
potentials as functions of the intermolecular separation. It is seen that while the spherically 
symmetric part is nearly identical for both potentials, the main anisotropic term, V022 = V202, 
responsible for the 00— >20 rotational excitation is smaller for the BMKP potential at small 
intermolecular separations. Our test calculations show that at energies lower than 0.1 eV, 
the discrepancy between the two results is mostly due to small differences in the coupling 
elements in the region of the van der Waals minimun, i.e., R > 3.0 au. The same also applies 
to the next higher order term, V224- On the other hand, Progrebnya and Clary^ found that 
the BMKP surface yields too high values for vibrational relaxation in H 2 (?; = 1)+H 2 (t> = 0) 
collisions. They attributed this to higher-order anisotropic terms in the BMKP PES that 
leads to preferential population of high rotational levels in the v ' — level after quenching 
(see Fig. 2 of the above reference). Calculations employing a modified version of the BMKP 
potential in which only the first two leading anisotropic terms of the interaction potential 
shown in Fig. 2 are retained gave results in better agreement with the experiment. 

The integral elastic cross sections of para-H2+para-H2 collisions as a function of collision 
energy is plotted in Fig. 3(a) for both the BMKP and DJ potentials. In the zero-energy limit, 
the elastic cross sections attain finite values in accordance with Wigner's law. The limiting 
value of the elastic cross section is 1.91 x 10~ 13 cm 2 and 1.74 x 10~ 13 cm 2 for the BMKP and 
DJ PESs, respectively. The comparable values of the limiting elastic cross sections on the 
two potentials is explained based on the nearly identical values for the spherically symmetric 
part of the interaction potentials for both surfaces (see Fig. 2). Both potentials exhibit a 
shape resonance at collision energies between 0.0002 and 0.0003 eV, arising from the / = 2 



partial wave. At higher energies, cross sections on the two potentials exhibit an oscillatory 
behaviour (see inset in the top panel of Fig. 3) which arises from interference between partial 
cross sections corresponding to different values of the total angular momentum quantum 
number J. Note that only even values of J are allowed for p-H 2 -p-H 2 collisions. Schaefer 
and Meyep^ have provided a detailed analysis of the oscillatory behavior of the elastic cross 
sections. 

In Fig. 3(b) we compare elastic cross sections from the present work on the DJ and BMKP 
surfaces with the theoretical results of Diep and Johnson 3 ^ and the experimental measure- 
ment of Bauer et al.— . It is seen that the results on the DJ potential give slightly better 
agreement with experimental data although the overall agreement between the experiment 
and theory is generally good. 

In Fig. 4 we compare cross sections from the present calculations on the BMKP surface 
for the 00 — > 20, 00 — > 22 and 00 — > 40 transitions with the wave packet results of Gatti et 
al. 33 and Lin and Guc 34 . The overall agreement between the CC and the wave-packet results 
is rather good, to within 10-15% for all cross sections, except for the 00 — > 22 rotational ex- 
citation for which the present results are 30-50% larger. The significant differences between 
the present results and the wave packet results for this transition is somewhat surprising 
considering the fact that the corresponding cross sections are about an order of magnitude 
larger than the 00 — > 40 transition for which we obtain good agreement with the wave packet 
results. We are confident that this is not a numerical error in our calculations as we have 
benchmarked our results against a new time-independent coupled channel code developed 
by Krems 4 ^ which reproduces our results up to several significant digits for all transitions 
shown in Fig. 4. The deviation of the CS wave-packet results of Lin and Guo from the 
present CC and full wave-packet data of Gatti et al. at higher energies may be attributed to 
Coriolis couplings. This evidently suggests that Coriolis coupling plays an important role at 
higher energies and that the rigid rotor approximation appears to hold well for this system 
for collision energies investigated in the present work. 

In Fig. 5(a), we compare rotational quenching cross sections for the 20 — > 00 transition 
evaluated using the BMKP and DJ potential with the 20 — > 20 elastic scattering cross 
section. Since the elastic scattering cross section on the two PESs is comparable, only the 
result on the DJ potential is shown. The resonance feature that occurs in all three cross 
sections just above 10~ 4 eV is due the I = 2 shape resonance discussed previously (See Fig. 



1). It is seen that the inelastic cross sections are much smaller than the elastic one at all 
energies shown in the Fig. 5(a). This suggests that it may be possible to cool rotationally 
excited H 2 molecules in the j = 2 rotational level by thermalizing collisions with ground 
state H 2 molecules (evaporative cooling) without significant trap loss, though the absence of 
an electric dipole moment makes it a difficult system to handle experimentally. For incident 
energies lower than 10 -5 eV, the quenching cross section varies inversely with the velocity in 
accordance with Wigner threshold behavior. As a consequence, the product of the relative 
velocity and the quenching cross section attains a finite value in the limit of zero incident 
kinetic energy, as illustrated in Fig. 5(b). The limiting value of the quenching rate coefficient 
is 2.4 x 1(T 13 cm 3 /s. 

In ultracold collisions where s-wave scattering dominates, elastic and inelastic scattering 
cross sections are conveniently expressed in terms of scattering lengths. The scattering length 
is real when only elastic scattering is present, but it becomes complex with the inclusion of 
inelastic channels^. The complex scattering length is given by aj 1 j 2 = ctj 1 j 2 — i(3j 1 j 2 where 
a and (3 are real and imaginary parts of the scattering length. The imaginary part of the 
scattering length /3 is related to the zero-temperature limit of the quenching rate coefficient: 
fejij,(T — ► 0) = 4^/3^2^///. For the DJ potential we obtain the values «oo = 5.88 A, 
«2o = 5.78 A, foo = 0.003 A, oi22 = 5.83 A, and $22 = 0.0023 A. The corresponding values 
for scattering on the BMKP potential are a o = 6.16 A, a 2 o = 6.16 A, (3 2 o = 0.00028 A, 
ct22 = 6.16 A, and $22 = 0.00071 A. It is seen that the real part of the scattering length 
remains practically unchanged for the three initial states indicating that no zero-energy 
resonances occur for any of the three initial states on either PES. The presence of zero-energy 
resonances (bound/quasibound states near channel thresholds) generally enhances the elastic 
scattering cross section at low energies. The smaller value of the inelastic quenching rates 
on the BMKP surface is also reflected in the values of /3 for all three initial states. 

In Fig. 6 we compare cross sections for 22 — > 00 and 22 — > 20 transitions obtained 
using the DJ potential with the CC calculations of Forrey^ based on the PES of Zarur and 
Rabitz 16 . The agreement is remarkably good considering that the DJ potential is derived 
from accurate ab initio calculations while that of Zarur and Rabitz is a model semi-empirical 
potential. The reproduction of the shape resonance near E = 10 -4 eV by the two potentials 
is a good indicator of the accuracy of the two potential surfaces. The solid curve in Fig. 6 
is the quenching cross section for 20 — > 00 transition on the DJ potential and it is seen that 



the cross sections for 22 — > 20 transition is larger at all energies. 

Further comparison between BMKP and the DJ potential is presented in Fig. 7 in which 
we provide energy dependence of the excitation cross sections for 00 — > 20, 22, 40, 42 and 
44 transitions. The general trend in all cases, except for 00 — > 40 at low energies, is that 
the BMKP potential yields smaller values of cross sections compared to the DJ potential. 
The differences get somewhat smeared out when the cross sections are integrated over a 
Boltzmann distribution of relative velocities of the two molecules to yield the rate constants. 
This is illustrated in Fig. 8 for the same transitions as given in Fig. 7. Comparison is also 
made to the rate coefficients computed by Flower— which show better agreement with the 
DJ results. 

IV. SUMMARY AND CONCLUSIONS 

We have performed quantum close-coupling calculations of elastic and inelastic rotational 
transitions in collisions of H 2 with H 2 using the two most recently published ab initio interac- 
tion potential energy surfaces (PESs) for the H4 system. The calculations span a wide range 
of energies (9-orders of magnitude) extending from the zero-temperature limit to about 2.0 
eV. Sensitivity of the results to details of the interaction potential is presented by computing 
real and imaginary parts of the scattering lengths for different initial rotational levels of the 
two colliding H 2 molecules. It is shown that the limiting elastic cross section is not very 
sensitive to the initial rotational levels of the two H 2 molecules although the inelastic cross 
sections strongly depend on the initial rotational level. 

We also showed that results obtained using the rigid rotor potential surface of Diep and 
Johnson^ are in close agreement with the experimental measurements of Mate et a/.— for 
00 — > 20 rotational excitation rate coefficient. However, the corresponding results obtained 
using the Boothroyd et al. 35 PES are a factor of ten smaller. This is quite significant 
because the Boothroyd et al. surface is generally believed to be the most accurate potential 
surface for the H4 system and it has been used in two most recent six-dimensional quantum 
mechanical calculations^*^ of rotational excitation in H 2 +H 2 collisions. The present study 
demonstrates that the BMKP surface will need to be reevaluated before it can be adopted 
in large-scale scattering calculations, especially for astrophysical applications. 
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FIG. 1: Rate coefficients for 00 — > 20 rotational excitation in H2+H2 collisions as a function of 
temperature. The solid and dashed curves are results obtained using the DJ and BMKP potentials, 
respectively. The solid squares are the experimental results and the dotted line with crosses is the 
theoretical calculation based on DJ PES by Mate et The results on the DJ potential agree 

very well with the experiment data whereas those obtained with BMKP do not. 
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FIG. 2: First three terms in the angular expansion of the intermolecular potential as functions of 
the intermolecular radial separation. The solid curves represent the DJ potential and the dashed 
curves denote the BMKP potential. 
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FIG. 3: Elastic cross section ctoo^oo as a function of collision energy. The solid and dashed curves 
represent the results on the DJ and BMKP surfaces, respectively. The dotted line is from Diep 
and Johnson 32 . The circles with error bars are measurements of Bauer et al4£. 




FIG. 4: Comparison of cross sections for 00 — > 20, 22 and 40 transitions. The calculations are based 
on the BMKP PES. The solid, dashed and dotted lines are the present close-coupling calculation, 
wave-packet calculation of Gatti et al—2, and Lin and Guo^i, respectively. 
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FIG. 5: Cross sections for 20 — > 00 rotational quenching as a function of collision energy, (a) 
Comparison with 20 — > 20 elastic cross section, (b) Energy dependent rate coefficients, i.e., relative 
velocity times cross section. 
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FIG. 6: Cross sections for 22 — > 00 and 20 quenching collisions as functions of collision energy 
plotted against 20 — > 00 cross sections. The calculations are based on DJ PES. The dotted lines 
are the results of Foney^. 




FIG. 7: Comparison of the cross sections for 00 — > 20, 22, 40, 42 and 44 transitions as functions of 
the collision energy. The solid and dashed lines represent DJ and BMKP PESs, respectively. 




FIG. 8: Similar to Fig. 7, except for the rate coefficient as a function of the temperature. 



